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We analyze the diffusion of a Brownian particle in a fluid under stationary flow. By using the 
scheme of non-equilibrium thermodynamics in phase space, we obtain the Fokker-Planck equation 
which is compared with others derived from kinetic theory and projector operator techniques. That 
equation exhibits violation of the fluctuation dissipation-theorem. By implementing the hydrody- 
namic regime described by the first moments of the non-equilibrium distribution, we find relaxation 
equations for the diffusion current and pressure tensor, allowing us to arrive at a complete description 
of the system in the inertial and diffusion regimes. The simplicity and generality of the method we 
propose, makes it applicable to more complex situations, often encountered in problems of soft con- 
densed matter, in which not only one but more degrees of freedom are coupled to a non-equilibrium 
bath. 



I. INTRODUCTION 



Liquid matter when subjected to the action of external forces or gradients exhibits peculiar characteristics which do 
not manifest in the absence of those external inputs. Its statistical mechanical properties present significant features 
as the appearance of long-range correlations and phase transitions, or the violation of the fluctuation-dissipation 
theorem (see for example Refs. [Q- 0). The aim of non-equilibrium statistical mechanical theories is precisely to 
explain the behavior of systems in such situations, away from equilibrium. 

These characteristics, far from being specific of pure substances as simple fluids Q, also manifest in complex fluids 
. Typical situations of transport in liquids or liquid-like systems involve the joint motion of liquid and solid or 
solid- like phases under the action of an external input. Many examples can be found in systems as polymers 
suspensions of neutral and field-responsive particles |ll]] and granular media under shear flow . The presence 
of shear flow significantly modifies transport properties, and may induce the appearance of new phases which otherwise 
would remain hidden (see for instance Refs. That is the reason why the influence of the shear flow in the 

dynamics is a subject which has received much attention in the last years, specially in the domain of soft condensed 
matter. 

Our purpose in this paper is precisely to analyze one of the simplest examples of a system coupled to a non- 
equilibrium bath, whose physical realization is a set of non-interacting Brownian particles moving in a fluid in station- 
ary flow. The suspended objects diffuse in and are convected by the fluid; their motion may eventually be influenced 
by the presence of external fields. These are the basic ingredients controlling the dynamics of the suspended phase. 
Since the Brownian objects are of mesoscopic nature, the dynamic description demands a mesoscopic treatment in 
terms of a probability distribution function. The evolution in time of this quantity is governed by Fokker-Planck and 
Smoluchowski equations. These equations constitute the basis for a mesoscopic description of the system and enable 
one to extract macroscopic information from the evolution equations for the moments of the distribution. It becomes 
then of primary importance to establish simple methods able to provide expressions of those equations in situations 
outside equilibrium. 

Kinetic equations of the Fokker-Planck type have been basically derived from kinetic theory of gases using the 
diffusion approximation in the Boltzmann equation ff|, |ll|[l7| , from the theory of stochastic processes through the 
master equation jl8| or by means of projection operator techniques [ p~9|j20| . It has also been shown that those kinetic 
equations can be derived from mesoscopic non-equilibrium thermodynamics (MNET) [pl| - |28f| . As in non-equilibrium 
thermodynamics, the basic point of this theory consists of assuming local equilibrium, which is performed at a more 
basic level: at mesoscopic level. This fact enables one to formulate a Gibbs equation in which the entropy, in 
accordance with the concept of Gibbs entropy, also depends on a 'density': the probability density. By applying the 
rules of non-equilibrium thermodynamics one obtains the entropy production and from it the corresponding linear 
laws between fluxes and forces. When using these laws in the continuity (balance) equation for the probability density 
one obtains the Fokker-Planck equation, pl] , p2|| . This is precisely the scheme we will adopt in this paper to analyze 
diffusion in stationary flow in both, the kinetic and the hydrodynamic regimes. 
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The paper is organized as follows. In Section II, we derive the Fokker-Planck equation describing Brownian motion 
under the presence of an arbitrary steady flow, in the framework of MNET. In Section III, we discuss the hydrodynamic 
description by constructing the hierarchy of moments of the distribution function, accounting for the macroscopic 
evolution of the system. These equations reveal the presence of inertial and diffusion regimes for the dynamics of the 
particles. In Section IV, we analyze the particular case of a shear flow. The corresponding Fokker-Planck equation 
is compared with the one obtained in |29| by means of kinetic theory. In the diffusion regime, we give explicit 
expressions for the pressure tensor and the viscosity. Finally, in the discussion, we summarize our main results and 
indicate potential applications of the formalism we have established to soft condensed matter systems. 



II. FOKKER-PLANCK DYNAMICS FROM MESOSCOPIC NON-EQUILIBRIUM THERMODYNAMICS 

We consider a dilute suspension of spherical particles of mass m, immersed in a liquid phase with constant density 
which acts as a heat bath. The whole system is subjected to conditions creating a stationary flow described by the 
velocity field 

v = v (r). (1) 

Since our main purpose is to analyze the effect of the velocity gradient on the dynamics of the particle, we will 
consider isothermal conditions neglecting 'viscous heating'. 

The mesoscopic nature of the suspended phase makes it necessary to analyze its dynamics by means of a Fokker- 
Planck equation accounting for the evolution of the distribution function, which may, in general, depend on the 
coordinates and momenta necessary to specify the state of the suspended objects. Since we assume no direct inter- 
actions between particles, the Brownian 'gas' will be described by means of the single-particle distribution function 
/(r, u, t), which depends explicitly on the position r , particle velocity u, and time t. 

The first step towards the obtention of the Fokker-Planck equation, is the formulation of the conservation laws for 
the 'gas' of suspended particles. In the absence of external body forces, the distribution function obeys the continuity 
equation 

§£ + V- Sf=~- 4 (2) 
at ou 

which introduces the current Js in phase space. The average of Eq.(||) with respect to the particle velocity u, leads 
to the macroscopic equation for the balance of mass, which can be written in the form 

Here p(r,t) is the density of the Brownian 'gas', given by 

p(r, t)=mj fdu, (4) 
v(r, t) is the average velocity of the Brownian particles defined through the expression 

pv(r, t) = m J ufdu. (5) 

and we have defined the total derivative as 

d d _ „ 

— = — +v-V. 6) 
dt dt K ' 

Our purpose is to obtain the equation governing the evolution of the probability density /, therefore we need to 
find out the explicit expression for the current J,-?. To this end, we will assume local equilibrium for which, entropy 
variations are given through the Gibbs equation [ pip2]| 

5s = ^Se + ^pSp- 1 - m J j;8c s du. (7) 

Here s(r, t) and e(r, t) are the entropy and total energy per unit of mass of the Brownian particles, respectively, p(f, t) 
is the hydrostatic pressure, p(u, r, t) is the non-equilibrium chemical potential per unit of mass, and = - is the 
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Brownian mass fraction pl[ . Notice that the term including the chemical potential in Eq.(0) is reminiscent of the 
corresponding one for a mixture in which the different components would be specified by the continuum 'index' u. 
Following the scheme of non-equilibrium thermodynamics |2l| , we will assume that Eq. (f7|) remains valid for changes 
in time and position into a mass clement followed along the center of gravity motion of the Brownian 'gas'. 

Since the particles undergo only translational motion, and in view of Eqs.(|^), (|^) and (Q), the remaining conservation 
law is the one for the energy. The presence of the external flow is responsible for the appearance in that equation of 
the term of 'viscous heating', giving rise to variations in the temperature field. To maintain isothermal conditions, and 
following previous ideas introduced in the implementation of the so-called 'homogeneous shear' [|o|J|T). we will assume 



the existence of a local heat source capable to remove the heat generated in the process. Under this assumption, the 
energy of the volume elements remains constant along their paths and its balance equation can be omitted in the 
subsequent analysis. 

The expression of the non-equilibrium chemical potential can be found through the Gibbs entropy postulate for the 
Brownian particles written as 



-k B J c s lnjf^du + s Leq . 



(8) 



Here k B is Boltzmann's constant and f Leq is the local equilibrium distribution function corresponding to the reference 
state described by the local Maxwellian with respect to the stationary convective flow jjj , 

f Leq (u,f 7 t) = e w^B-i(S-i?o) 2 ], ( 9 ) 

In this expression /is(r, t) is the local equilibrium chemical potential of the Brownian particles. The local equilibrium 
entropy for the Brownian particles per unit of mass is given by 

s Leq = ~ J csVBdu+^e. (10) 

After applying the total derivative in Eqs. (0), (||) and ([h]), and comparing the corresponding expressions for the 
time variations of the entropy, we arrive at the desired expression for the non-equilibrium chemical potential 

H= —lnf+\(u-v ) 2 . (11) 

This expression, together with the equation for the Brownian mass fraction 

^ = - p V.v -V.f(u-v )-JL.J u> (12) 
dt ou 

obtained by combining the continuity equation (Q) and the balance of mass (|3|), is now substituted into Eq.(^|). The 
resulting expression is integrated by parts over the velocity space assuming that the fluxes vanish at the boundaries. 
One then obtains the entropy balance equation in the form 

ds - , 

p— = -V • J s + a, (13) 

where the entropy flux is given by 

J s = -k B (u-v )f(lnf-l)dti- — (u - v )f(u - v ) 2 du, (14) 



and the entropy production is 



a = -— / J u ■ T-du - — / J • (u - v ) ■ Vvadu. (15) 



T J du T 

This quantity consists of two contributions of the type flux-force pair: the first one arises from the diffusion process in 
M-space, whereas the second is due to the presence of the convective flow J = (u — vo)f. Since both contributions are 
vectors, the fluxes couple to the two forces giving rise to cross-effects. Following the non-equilibrium thermodynamics 
rules, we can establish linear phenomenological relationships between fluxes and thermodynamic forces. Assuming 
locality in u-space, for which only fluxes and forces with the same value of u are coupled, the expressions for the 
currents are the following 
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■Az = --^Lm ■ tt^ - T^Litf [u - v ) ■ \/v , 
1 Ou 1 



(16) 



J = -7^ L rr ■ [U - Vq) ■ VVo - — LfS ' 

1 T ou 



(17) 



where Las, ^™ and £?r Onsager coefficients. These coefficients may in general depend on the imposed 

velocity gradient. In this case they satisfy generalized Onsager relations |32] ]. in which time-reversal symmetry must 

also be applied to the external driving. In the case we are considering, we have Lap — —L^. Defining now the tensors 



and using the expression for the chemical potential ([ll]), the fluxes can be recast in the form 

Ju = ~{u - v ) ■ [a + Vv ■ e\f a ■ — , 

to Ou 



(18) 



(19) 



J 



-(u-vo) 



Vv • £ - e 



f + — e ' u^- 

m cm 



(20) 



By substituting Eq. ( |19|) into the continuity equation (g) for the single-particle distribution function, we finally obtain 

(21) 



(it - w ) • a + Vwo • H a ■ — 

m ou 



which constitutes the Fokker-Planck equation describing the evolution of the non-equilibrium single-particle distribu- 
tion function. The fact that the coefficients appearing in the equation are tensors, reflects the anisotropy of the system 
induced by the imposed external flow. Moreover, this equation exhibits the fact that the flow breaks the Einstein 
relation by adding a term which depends on the imposed velocity gradient. This breaking constitutes a proof that 
the fluctuation-dissipation theorem cannot be applied when fluctuations take place around the steady state. 

The Fokker-Planck equation we have obtained can be compared with the ones derived by means of different methods. 
In Ref. [p9l , authors found a similar Fokker-Planck equation for the particular case of a shear flow. In Section IV, we 
will discuss that similarity in more detail. Following time-dependent projector operator techniques, a Fokker-Planck 
equation similar to ours was derived in |33| for the translational modes of a Brownian particles moving in a flowing 
bath under a temperature gradient. Since diffusion in a bath under temperature gradient was studied previously in 
|p2| in the framework of MNET, the relevant point to emphasize here is the fact that the term including the velocity 
gradient enters the Fokker-Planck equation as an external force, in a similar way as in the equation obtained in p3| . 



III. THE HYDRODYNAMIC EQUATIONS 

In the absence of direct interactions, the single-particle distribution function provides the complete description of 
the system at mesoscopic level. Macroscopically, the description must be carried out in terms of the moments of the 
distribution function, which are related to the hydrodynamic fields: density, momentum and pressure tensor. 

The density defined through Eq. (|j) corresponds to the zero-order moment. The first-order moment has been 
defined in (||). The second moment centered about the average velocity v of the Brownian 'gas' corresponds to the 
pressure tensor 



P = to J (u-v)(u- v)fdu. (22) 
Finally, the third centered moment, which is related to the flux of kinetic energy and stress, is given by 



Q = m (u — v)(u — v)(u — v)fdu. 



(23) 
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The set of evolution equations for these moments can be obtained by using the Fokker-Planck equation (^T|) in the 
definitions of the conserved quantities, after performing the time derivative and the required integrations in velocity 
space pH]. The evolution equations for the three first moments are, respectively, the continuity equation 

| = "V-M (24) 



the balance of momentum 



dv ~—* ~—* 

p- n +V-P=-p(v-v Q )-C, (25) 



dt 

and the equation for the evolution in time of the pressure tensor 

j t P + 2(P ■ Vv) s + PV-v + 2(P- C) s = 

= ^ p( T-v4 (26) 
m 

where C = [a + Vt/o • e] and an upper s means symmetric part of a tensor. In a similar way, we could derive 
the evolution equations for the higher-order moments of the distribution, which constitute a coupled hierarchy of 
hydrodynamic equations p5] , |34[ . 

Notice that in Eq. (|2^) the right hand side term can be identified with the hydrodynamic force exerted by the fluid 

on the particle, with C playing the role of the friction constant. That friction constant establishes the characteristic 
relaxation time scale for the velocity and it is usually very large. For instance, for a mesoscopic particle of radius 
a ~ 10~ 5 cm, moving in a quiescent liquid of viscosity rj ~ 10 _2 P, the Stokes formula can be applied giving for the 
friction constant (3 = ^112. ^ io 8 s _1 . Consequently, the discussion of the behavior of the system may be carried 

out by expanding the hierarchy of evolution equations for the moments in powers of C . From the hierarchy of 

moments, one can easily realize that the ith moment introduces corrections which are of order C for even 

-* — i 

moments, and C for the odd ones. If we drop out the terms arising from the third and higher moments, we are 

^-3 3-3 

neglecting corrections to the diffusion equation of order C . Notice that in order to retain C corrections, one has 
to include not only the third-order moment equation, but also the equation for the forth moment which is of the same 
order. 

The set of equations (|2^)-([2^) then govern the hydrodynamic behavior of the Brownian 'gas' immersed in a fluid 

moving with velocity profile vq. The elements of the tensor C constitute characteristic time scales whose existence 
motivates the separation of the dynamics into two well-differentiated regimes: an inertial regime for t <C {C 
characterized for the relaxation of the variables towards the diffusion regime, which is achieved for t 3> (C -1 )^ . 
Both regimes will be discussed in the following subsections. 



A. Inertial regime 

In order to discuss the inertial regime, it is convenient to rewrite Eq. (^6|) for the evolution of the pressure tensor 
in the following way 

where we have defined the matrix of relaxation times 

?i = + W+ (V-vjf]- 1 , 

3-3 

and we have neglected corrections to the diffusion equation of orders C and higher. In an analogous way, from the 
evolution equation for the momentum, Eq. (p5|) , we obtain 

^jT + Jd ■ ?^ - Pvo ■ 3 - V • P, (28) 
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where Jd = pv and 

7 2 = 0+i(V -v)}- 1 , 

is the corresponding matrix of relaxation times. 

The previous equations describe the inertial regime in the dynamics of the Brownian particles subjected to stationary 
flow. This regime holds for times small enough compared to the characteristic relaxation times, identified with the 
components of the matrices t\ and T2 ■ 

Our description of the inertial regime is consistent with the generalized hydrodynamic description in which the 
diffusion coefficient depends on the wave vector |35| . This property has been shown in ]2q] for the case of a quiescent 
liquid. 



B. Diffusion regime. 

For times larger than any characteristic relaxation time, t 3> {C the system enters the diffusion regime. In 
such a regime, the dynamics becomes well described by a Smoluchowski equation for the density distribution p in the 
configuration space. 

In order to find the Smoluchowski equation, we will first discuss the diffusion approximation in the evolution 
equations for the pressure tensor ( p7| ) and the momentum (^8|). Both equations involve inertial, friction and velocity 
gradient time scales. For t 3> (C~^lj, time derivatives can be neglected when compared with terms proportional to 

C. Notice also that the term V • v is essentially a time derivative, as follows from Eq. (||), and can accordingly be 
neglected. Taking this considerations into account, the equations for the pressure tensor and the momentum then 
reduce to 

(P ■ dy + (P • w) s = — P a , (29) 

m 

and 

Jd =pvo- (v-pyd , (30) 

respectively. The first equation represents a linear set of coupled algebraic equations for the components of the 

pressure tensor P in terms of the components of the tensors C, V?7 and a. Once the explicit expression for the 
pressure tensor has been found, it must be introduced into the momentum equation ( |30| ) yielding the constitutive 
equation for the diffusion current Jd- The Smoluchowski equation is then obtained after substituting Jd into the 
balance of mass (pi|). 

In the following Section we will apply the previous scheme to the particular case of a shear flow. 



IV. DIFFUSION IN A SHEAR FLOW 



Inherent to non-equilibrium thermodynamics is the fact that it cannot provide explicit expressions for the phe- 
nomenological (transport) coefficients which must be borrowed from other theories. In our case, the still unspecified 
quantities are a and e. When a specific velocity profile is given, expressions for those tensors can be obtained from 
kinetic theory or hydrodynamics. 

The formalism developed in the previous sections is valid for an arbitrary stationary velocity field. In this section 
we will focus our discussion on the particular case of Brownian motion in a shear flow. Our first task will be the 
identification of the phenomenological coefficients. As pointed out before, the second term on the right hand side 
of equation (|2^) can be identified with the force per unit of mass exerted on a suspended particle by the host fluid, 

where C plays the role of the friction tensor. In general, such a density force can be calculated from hydrodynamics, 
(see for instance Refs. |}6|,[37]]). For the shear flow case, the friction coefficient contains in general linear and quadratic 

contributions in (a, where £ is the inverse penetration length of the perturbations £ = ( J) 2 , with 7 being the shear 
rate and v the kinematic viscosity. For a Brownian particle under moderate shear flow the term £a is very small, 

consequently the friction tensor can be approximated by C ~ 01. This identification leads to the following expression 
for the tensor a, 



G 



a = (31 — S7vq ■ e. 

By introducing the previous expressions into the Fokker-Planck equation ( pl| ) , we finally obtain 

d f__ k R T 



at 



u-Vf= — -\(3{u-v )f 



<BJ- i df_ 
m du 



(31) 



(32) 



As concluded also from Eq.fl2l|), the particular form of Eq.(32) implies that the fluctuation-dissipation theorem is no 
longer valid when the fluid (heat bath) is sheared. Notice that the term which invalidates that theorem is proportional 
to the velocity gradient, or to the inverse penetration length squared. The theorem remains applicable in the case in 
which the friction coefficient contains a correction proportional to £ 38 1 . 

The Fokker-Planck equation (^2|) is similar to the corresponding one obtained in Ref. |^9| from kinetic theory. 
By expanding the collision operator in the mass ratio between the fluid and Brownian particles, authors derived a 
Fokker-Planck equation in which the diffusion tensor contains a correction to the Einstein formula proportional to 
the stress tensor of the fluid. Since this quantity is proportional to the velocity gradient we conclude that the form of 
the diffusion tensor they find is similar to our expression (|3l| ) , which allow the identification of the tensor e. 

With the explicit form of the Fokker-Planck equation (|32|) in mind, our purpose is now to discuss the macroscopic 
evolution of the system. Following the procedure indicated in subsection III.B, we find that for a shear flow, in the 
diffusion regime, the expression for the Brownian pressure tensor is 



P 



knT 



-p 



1 



(33) 



From this equation, we can conclude that Brownian motion of the particles contributes to the total pressure tensor 
of the suspension in two forms. The first contribution is the well-known scalar kinetic pressure given by 



k B T 



P = 



-P< 



(34) 



which is the equation of state for the ideal Brownian 'gas'. The second contribution comes from the irreversible part 
n of the Brownian pressure tensor, which can be written in the form 



n : 



D p (i+e)-Vv 



(35) 



where Dq = is the diffusion coefficient of a particle when the liquid is at rest. This last equation defines the 
Brownian viscosity tensor 



Vb = D p(l + e), 



(36) 



which contains the 'Brownian' viscosity D p |39|] , and the contribution due to the coupling with the non-equilibrium 
bath which is proportional to <f. Eqs. ( |35| ) and d36| ) accounts for the contributions to the irreversible part of the 
pressure tensor and shear viscosity coefficient of the suspension. No contribution to the bulk viscosity has been found 
since the shear flow is incompressible. 

Following the steps described previously, we can obtain the diffusion current Jd in the limit t 3> {C~ l )ij, in which 
Eq.@ has the form 



Jd = pv - /T *V ■ P- 

Introducing the expression for the pressure tensor (133) into the last equation we obtain, 



Jd = pvo - D ■ Vp - pV • D. 
Here we have defined the diffusion tensor as 



D = Dn 



1 



/r 1 ((i + f)-w o y 



(37) 



(38) 



(39) 



An important consequence of this result is that the presence of a shear flow modifies the diffusion current with respect 
to the case of a quiescent liquid. Substituting Eq.(|38|) into the balance of mass, Eq.(p4]), we finally obtain the diffusion 
equation for the Brownian particle 
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^ = -V • pv + V • • Vp) + V • (pV • D). (40) 

This equation is a generalization of the usual Smoluchowski equation, in the sense that it contains the anisotropic 
diffusion coefficient (pa) which depends on the imposed velocity gradient. For the particular case of the shear flow, the 

velocity gradient is a constant, therefore D does not exhibit spatial dependence and we recover the usual Smoluchowski 
equation 

% = -V • pv + V • ■ Vp). (41) 

This equation coincides formally with the one found in Ref. |fi"o|j . The difference between both equations lies in the 
form of the diffusion tensor. In the cited work, author starts from a Fokker-Planck equation which satisfies the 
fluctuation-dissipation theorem. In the diffusion regime author finds, consistently with the approximation, that the 
diffusion tensor and the mobility of the Brownian particles in the fluid are related each other through Einstein's 
formula. 



V. DISCUSSION 



In this paper we have analyzed the dynamics of a suspension of Brownian particles in a non-equilibrium situation 
resulting from the action of an externally imposed velocity gradient. We have applied the method of mesoscopic 
non-equilibrium thermodynamics to study the dissipation in phase space related to the underlying diffusion process 
of the probability density of the particles. 

In MNET, local equilibrium is assumed at mesoscopic level. A Gibbs equation is then proposed in which the entropy 
depends on the probability density: the pertinent density in phase space, according to the concept of Gibbs entropy. 
By applying the rules of non-equilibrium thermodynamics, one obtains the entropy production of the system, which 
enables one to derive the expression for the diffusion current in phase space and, consequently, to obtain the Fokker- 
Planck equation for the single-particle distribution function. This expression exhibits violation of the fluctuation- 
dissipation theorem whose origin is precisely the presence of the external gradient. This feature, commonly found in 
the wide class of driven-diffusion systems fjJI^], has also been reported in slow relaxation process of glassy systems 

The hydrodynamic level of description is accomplished from the evolution equations for the first moments of the 
distribution function, which can be obtained through the Fokker-Planck equation. The time evolution of the moments 
include relaxation equations for the diffusion current and the pressure tensor, whose form permits to elucidate the 
existence of inertial (short-time) and diffusion (long-time) regimes. In the diffusion regime, the mesoscopic description 
is carried out by means of a Smoluchowski equation. In this regime, the equations for the moments coincide with the 
differential equations of non-equilibrium thermodynamics. The equations for the moments obtained from the Fokker- 
Planck equation in the framework of MNET extend the domain of applicability of non-equilibrium thermodynamics 
to shorter time scales. These equations can be reformulated in terms of transport coefficients which depends on the 
wave vector, in accordance with generalized hydrodynamics p5j. 

Our results can be compared with others obtained previously using different theories. The Fokker-Planck equation 
we have derived is similar to the one proposed in |29| for the case of a shear flow. As in our case, these authors 
show that when the system is sheared, the diffusion coefficient is no longer given by Stokes-Einstein law. It contains 
a correction proportional to the pressure tensor which is basically of the same nature as the one we obtain, i.e., 
proportional to the velocity gradient. This fact clearly shows violation of the fluctuation-dissipation theorem due to 
the presence of the external input necessary to maintain the stationary state. Similar conclusions are obtained in 
p3| , where the general case in which a stationary flow and a temperature gradients act simultaneously is studied by 
means of time-dependent projector operator techniques. 

Having discussed the comparison of our results with others coming from statistical mechanical theories, it remains 
to analyze them in the framework of thermodynamical theories dealing with systems outside equilibrium, in particular 
with extended irreversible thermodynamics. This theory provides hydrodynamic equations which also contain relax- 
ation terms for the pressure tensor and the diffusion current. Along our analysis of this and other cases of systems 
outside equilibrium (see Refs. |^4|,^5],^7 42 1), we have shown that by simply using non-equilibrium thermodynamics 
or its extension to the mesoscopic domain (MNET) , we are able to completely characterize the evolution of systems 
outside equilibrium. This fact question the need of using generalized entropies p3[ | depending on non-thermodynamic 
variables: the fluxes, which constitutes the cornerstone of extended irreversible thermodynamics. The application 
of the well-established non-equilibrium thermodynamics postulates as indicated in pl[ | , suffices to provide a general 
scheme under which non-equilibrium processes of macroscopic and mesoscopic nature can be treated. 
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Far from being specific of the case of a suspension of non-interacting Brownian particles, we have treated in this 
paper, the method we have presented could systematically be applied to analyze the dynamics of soft condensed 
matter systems under shear flow. For example, the case in which direct and hydrodynamic interactions among 
particles become important can be worked out along the lines indicated in pi) , for the case of a quiescent liquid. If 
the suspended objects are deformable or need of additional parameters to characterize their state, as may occur in 
polymers |9) or liquid crystals the set of variables of the distribution function has to be enlarged. One would 

obtain the corresponding Fokker-Planck equation and from it the evolution equations for the moments which define 
the hydrodynamic or generalized hydrodynamic regimes. In all these cases, kinetic and hydrodynamic equations can 
be derived following the method of mesoscopic non-equilibrium thermodynamics. 
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